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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  65807F. 
The  results  of  the  research  were  obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup 
Corporation  Company),  operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force 
Station,  Tennessee,  under  ARO  Project  Number  P32A-K4A.  Analysis  of  the  data  was 
completed  August  31,  1977,  and  the  manuscript  was  submitted  for  publication  on 
December  14,  1977. 

The  author  wishes  to  acknowledge  the  contribution  of  R.  W.  Butler  and  W.  C 
Armstrong,  ARO,  Inc.,  who  supported  and  aided  the  research  effort. 
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1.0  INTRODUCTION 

The  developed  spin  is  classically  defined  as  an  equilibrium  motion  in  which  an 
airplane  is  descending  toward  earth  at  some  angle  of  attack  between  stall  and  90  deg 
while  rotating  about  the  vertical  or  near-vertical  axis.  The  nature  of  this  developed  spin 
can  be  steady  or  oscillatory.  From  a pure  definition  sense,  the  steady  spin  occurs  when 
the  aerodynamic,  inertial,  and  gyroscopic  moments  of  the  airplane  are  balanced.  An 
oscillatory  spin  occurs  when  the  moment  equations  of  motion  are  not  balanced  and  thus 
are  simultaneously  causing  repetitive  motions  in  both  frequency  and  amplitude.  The 
airplane  flat  spin  is  characterized  as  a steady  spin  and  is  accompanied  by  extreme  angles 
of  attack,  high  rotation  rates,  and  small  spin  radius.  In  reality,  those  flight-experienced 
spins  that  are  characterized  as  steady,  flat  spins  are  generally  accompanied  by  a 
superposition  of  small  oscillations  in  all  rates  and  angles  on  the  steady  motions. 

Flight  experience  has  shown  that  the  developed  flat  spin  can  be  entered 
inadventently  in  high  angle-of-attack  maneuvers  and  that  recovery  can  be  difficult,  if  not 
impossible.  For  these  reasons,  much  effort  has  been  expended  to  better  understand  and 
predict  airplane  spin  behavior. 

Successful  analytical  simulation  of  airplane  spinning  motion  has,  in  most  cases,  been 
achieved  only  after  flight  test  data  were  acquired.  With  the  aid  of  flight  data,  the 
aerodynamic  data  matrix  is  altered  to  provide  good  flight/simulator  spin  correlation.  This 
after-the'-faet  adjustment  is  required  because  the  true  aerodynamics  of  the  airplane 
spinning  motion  cannot  be  defined  properly  in  conventional  static  wind  tunnel  testing. 
Large  angular  motions  and  separated  flow  conditions  combine  to  produce  non  linearities 
in  the  aerodynamic  forces  and  moments  acting  on  the  airplane  in  the  spin.  The  severity 
of  these  nonlinearities  is  dependent  on  the  configuration  of  the  airplane  and  is  the  prime 
factor  in  determining  the  quantity  of  aerodynamic  data  required  for  accurate  simulation 
of  spinning  motions. 

Attempts  to  analytically  model  departure  and  spin  entry  characteristics  using 
measured  static  and  forced-oscillation  aerodynamic  data  (conventional  data  matrix)  were 
undertaken  by  the  National  Aeronautics  and  Space  Administration  (NASA)  (Refs.  1 and 
2).  The  NASA  investigation,  using  the  six-degree-of-freedom  equations  of  motion,  did  not 
result  in  a good  representation  of  the  departure  and  spin  entry  flight  motions.  The  results 
did  show  that  yawing  moments  produced  at  high  angles  of  attack  were  nonlinear  with 
respect  to  the  rotation  rate  of  the  airplane.  Simulated  flight  motions  using  measured  large 
angle  static  aerodynamics  did  produce  more  representative  motions  in  the  stall  region 
(Ref.  3). 
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Because  of  the  apparent  inability  to  realistically  simulate  spin  motions  using  the 
conventional  aerodynamic  data,  NASA  Langley  began  a program  to  measure  aerodynamic 
data  on  airplane  models  under  constant  rotating  conditions.  This  program  resulted  in  the 
development  of  a rotation-balance  apparatus.  The  resulting  data,  measured  at  steady 
rotating  conditions,  were  used  initially  for  spin  tunnel  test  correlations  (Ref.  4). 
Additional  data  obtained  from  a similar  rotation-balance  apparatus  (Ref.  5)  were  used  in 
an  analytical  study  of  aircraft-developed  spin  and  spin  recovery  characteristics  (Ref.  6). 
The  configurational  dependence  of  the  nonlinearities  in  the  aerodynamics,  with  respect  to 
rotation  rate,  was  demonstrated  by  the  data  presented  in  Refs.  5 and  6. 

A refined  method  for  obtaining  more  accurate  rotation-balance  aerodynamic  data, 
employing  larger  models,  was  developed  by  NASA.  The  data  obtained  from  this  new 
technique  have  been  used  to  analytically  simulate  representative  spinning  motions  of  a 
1/ 10-scale  radio-controlled  airplane  model  (Refs.  7 and  8). 

The  objective  of  the  work  presented  in  this  report  is  to  investigate  the  usefulness  of 
rotation-balance  wind  tunnel  data  in  computer  simulation  of  full-scale  airplane  spinning 
motion.  The  investigation  uses  a nonlinear  six-degree-of-freedom  digital  computer  program 
with  aerodynamic,  inertial,  and  thrust  data  for  the  F-4  series  fighter  aircraft.  Simulated 
spins  are  generated  utilizing  (1)  static  and  forced-oscillation  data,  and  (2)  static, 
forced-oscillation,  and  rotation-balance  data.  These  simulated  spins  are  in  turn  compared 
with  time  histories  of  F-4E  flight  test  spins.  Specifically,  the  discussion  and  evaluation  of 
the  application  of  rotation-balance  data  in  airplane  motion  simulation  will  be  limited  to 
the  steady,  flat  spin. 


2.0  ANALYSIS 


2.1  DESCRIPTION  OF  AIRPLANE 

The  selection  of  an  airplane  configuration  to  be  used  in  this  investigation  was 
influenced  by  several  factors.  A current  fighter  airplane  configuration  for  which  high 
angle-of-attack  wind  tunnel  data  were  available  was  desired.  Specifically,  a comprehensive 
aerodynamic  data  matrix  including  static,  forced-oscillation,  and  rotation-balance  data  was 
deemed  necessary.  In  addition,  flight  test  documentation  of  the  basic  spin  modes  and 
availability  of  these  data  were  required.  Recent  wind  tunnel  tests  conducted  at  Langley 
Research  Center  produced  a large  amount  of  measured  high  angle-of-attack  aerodynamic 
data  on  the  F-4  configuration  which  was  documented  in  Refs.  9 and  10.  The  Air  Force  has 
also  conducted  extensive  full-scale  flight  tests  investigating  the  stall,  departure,  and  spin 
characteristics  of  the  F-4E  configuration.  On  the  basis  of  the  vast  amount  of  aerodynamic 
data  and  flight  test  results  available,  the  F-4  configuration  was  selected  for  this  investigation. 
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■A  three-view  sketch  of  the  general  F-4  series  airplane  configuration  is  shown  in  Fig. 
1.  The  mass;  inertial,  and  geometric  characteristics  used  in  this  investigation  are  listed  in 
Table  1.  These  characteristics  are  for  a clean  airplane  with  no  external  stores  other  than 
.sparrow,  missiles  at  weapon  stations  4 and  6.  The  F-4E  characteristics  were  chosen  in 
order,  to  be  consistent  with  the  flight  test. 

2.2  AERODYNAMIC  DATA 

' The  major  part  of  any  analytical  investigation  which  calculates  the  steady-state  spin 
motion  of  an  airplane  lies  in  defining  the  aerodynamic  model  to  be  used.  The  present 
study  utilizes  three  types  of  measured  force  and  moment  aerodynamic  data:  conventional 
static,  .forced-oscillation  dynamic,  and  rotation-balance. 

All  aerodynamic  data  used  in  this  investigation  were  compiled  from  tests  conducted 
in  the  Langley  full-scale  tunnel  at  a speed  of  approximately  100  ft/sec  (Mach  No.  *=  0.1), 
which  corresponds  to  a Reynolds  number  of  about  1.35  x 106  based  on  wing  mean 
aerodynamic  chord.  Linear  interpolation  was  used  on  all  aerodynamic  table  data,  unless 
otherwise  noted. 

Conventional  Data  Matrix  - The  usual  method  for  representing  the  aerodynamic 
characteristics  of  an  airplane  is  a combination  of  static  and  forced-oscillation  wind  tunnel 
data.  The  static  data  used  were  obtained  from  Ref.  9 and  included  the  following 
coefficients: 

Cx  (o,  /3,  SH)  ; Cj,  (a,  j8,  SH)  ; Cm  (a,  1 3,  51{) 

(a  , £ , 8.)  ; Cg  (a,  /3 , Sa)  ; Cn  (a  , /3 , 8 a) 

Cy  (a,  /3 , 8)  ; Cg  (a,  jS,  5r)  ; Cn  (a,  (3,  8f) 


..  These  coefficients  were  measured  at  angles  of  attack  from  -10  to  110  deg  in 
increments  of  5 deg  and  sideslip  angles  of  0,  ±5,  ±10,  ±15,  ±20,  ±30,  and  ±40  deg. 

■ . i 

The  forced-oscillation  derivatives  were  obtained  from  Ref.  10  and  included  the 

t 

following  damping  parameters: 


Cm  ~ Cm-;  Cv  * C v A sin«i  Cv  - CvAC0Sa 


n ■ v yo 

a ' p J p 


' j8 


C „ + C . ; C„  + C . sin  a ; C„  - C . cos  a ; 

z„  za  nP  /3  "r  n|8 

G x q + CV  C«P  + CZpsinai  C£r  “ 
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It  should  be  noted  that  the  measured  coefficients  contain  terms  in  a and  j3  which 
arise  from  the  forced-oscillation  measuring  technique.  For  this  investigation,  these  terms 
are  assumed  to  be  zero,  as  is  conventionally  done,  and  the  measured  damping  terms  are 
used  as  pure  derivatives  in  p,  q,  and  r rates. 

These  forced-oscillation  coefficients  were  measured  at  frequencies  of  0.7,  1.0,  and 
1.3  cycles  per  second  and  amplitudes  of  ±5  and  ±10  deg  over  an  angle-of-attack  range 
from  -10  to  110  deg  at  5-deg  increments.  Also,  the  "r"  derivatives  were  measured  with 
elevator  deflections  of  0 and  -21  deg  for  the  three  frequencies  at  an  amplitude  ratio  of 
±5  deg  over  the  angle-of-attack  range.  Only  the  data  corresponding  to  the  reduced 
frequency  parameter,  K(£2  b/2V  or  SI  C/2V),  obtained  from  an  analysis  of  the  F-4  spins 
(Ref.  11)  and  from  the  static  wind  tunnel  data  (Ref.  9)  using  a first-order 
approximation  [(cosp  = f(Cm  a /ly ),  to  = f(Cnj}/Iz)]  were  incorporated  into  the  simulation 
data  matrix.  The  longitudinal  and  lateral/directional  damping  characteristics  incorporated 
were  for  a frequency  of  0.7  cycles  per  second  and  an  amplitude  of  ±5  deg. 

Rotation-Balance  Data  Matrix  - The  rotation-balance  wind  tunnel  data  were  obtained 
with  the  wind  tunnel  velocity  vector  coincident  with  the  airplane  model  rotation  vector.  A 
sketch  of  the  rotation-balance  test  apparatus  used  in  the  full-scale  wind  tunnel  at  NASA 
Langley  is  shown  in  Fig.  2.  The  following  coefficients  were  measured  with  the 
rotation-balance  apparatus:  Cx,  CZ)  Cm  , Cg,  Cn,  and  Cy. 

These  total  coefficients  were  measured  at  nondimensional  rotation  rates  (Sib /2V)  of  0, 
±0.049,  ±0.098,  ±0.196,  ±0.245,  and  ±0.294,  pitch  angles  of  55  to  90  deg  in  5-deg 
increments,  and  bank  angles  of  0,  5,  and  10  deg.  The  tunnel  roll  and  pitch  angles  are 
related  to  the  airplane  angle  of  attack  and  sideslip  by 

a - tan-'  (tan  6 cos  o) 

/3  = sin-'  (sin  0 sin  <f>) 

The  incremental  rotation-balance  data  used  in  this  investigation  included  data 
measured  with  all  controls  neutral,  elevator  control  at  -21  deg,  and  with  left  and  right 
pro-spin  controls  deflected  (left:  8h  = -21  deg,  5a  = -30  deg,  5r  = 30  deg;  right:  5h  = -21 
deg,  5a  = 30  deg,  5r  = -30  deg).  To  represent  the  aerodynamic  damping  effect  of  oscillations 
superimposed  on  a steady  rotation,  forced-oscillation  derivatives  were  used  with  the 
rotation-balance  data. 
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Significant  magnitudes  of  rotation  rate  are  known  to  develop  in  fighter  airplane 
stall/post-stall  gyrations  below  55  deg  angle  of  attack.  The  incremental  rotation-balance 
data  effects  at  55  deg  angle  of  attack  were  extrapolated  to  30  deg  angle  of  attack  for 
possible  excursions  in  this  region  due  to  control  perturbations. 

2.3  METHOD  OF  ANALYSIS 

Airplane  steady-state  spins  are  generally  characterized  by  sustained  turning  motion  at 
a stalled  angle  of  attack.  To  adequately  represent  these  motions  in  an  analytical 
investigation  requires  the  complete  six-degree-of-freedom  equations  of  motion.  A 
comprehensive  matrix  of  measured  nonlinear  aerodynamic  data  are  also  required  for  these 
high  angle-of-attack  flight  regions.  The  following  discussion  describes  the  equations  of 
motion  and  the  aerodynamic  data  models  used  in  the  investigation. 

The  differential  equations  of  motion  are  presented  in  Appendix  A (see  A.1).  Also 
presented  are  equations  defining  parameters  required  by  the  six-degree-of-freedom 
equations.  The  method  by  which  the  aerodynamic  force  and  moment  data  were 
combined  to  evaluate  the  inclusion  of  rotation-balance  data  in  spin  simulation  will  now 
be  described.  (Further  discussion  may  be  found  in  Appendix  B.)  Three  aerodynamic 
models  were  formulated  from  the  previously  described  force  and  moment  data.  Aero 
Model  1 represents  the  conventional  technique  which  does  not  include  the  aerodynamic 
forces  and  moments  acting  on  the  spinning  airplane  because  of  steady  rotational  flow  and 
does  not  restrict  the  dynamic  oscillatory  derivatives  to  the  oscillatory  component  of  the 
total  angular  rates.  Aero  Model  2 represents  the  conventional  technique  in  which  the 
dynamic  oscillatory  derivatives  are  restricted  to  the  oscillatory  component  of  the  total 
angular  rates.  Aero  Model  3 represents  the  aerodynamic  model  based  on  the 
rotation-balance  data  in  which  the  dynamic  oscillatory  and  rotation-balance  data  are 
restricted  to  the  oscillatory  and  steady  rotation  components  of  the  total  angular  rates, 
respectively. 

The  initial  conditions  and  control  deflections  used  in  the  analytical  simulations  were 
obtained  from  Ref.  1 1 and  are  documented  in  Table  2.  These  initial  conditions,  along 
with  the  mass  and  inertia  characteristics  and  each  of  the  aerodynamic  models,  were  input 
into  the  computer  program.  The  program  was  numerically  intergrated  with  time  by  using 
a fourth-order  Runge-Kutta  integration  algorithm  with  a fixed  integration  stepsize  (Ref. 
12).  An  examination  of  the  computed  motions  for  each  aerodynamic  model  then 
indicated  the  success  of  the  simulation  attempt. 
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3.0  RESULTS  AND  DISCUSSION 

The  simulation  of  a steady,  flat  spin  was  performed  using  the  three  aerodynamic 
data  models  described  in  Section  2.3.  A flight  test  spin  characterized  as  a flat  spin  was 
selected  as  a correlation  base  because  of  the  proximity  of  the  velocity  and  spin  rate 
vectors.  The  relationship  of  the  velocity  and  spin  vectors  of  the  chosen  spin  more  closely 
agreed  with  the  conditions  under  which  the  rotation-balance  wind  tunnel  data  were 
obtained.  This  particular  spin,  although  characterized  as  a flat  spin,  was  accompanied  by 
small  oscillations  in  all  rates  and  angles.  These  oscillations  were  repetitive,  with 
amplitudes  of  less  than  ±20.0.  The  initial  conditions  used  in  the  simulation  were 
obtained  from  a time  history  presented  in  Ref.  1 1 documenting  the  F-4E  stall/near-stall 
characteristics.  The  necessary  parameters  required  for  the  simulation  were  extracted  from 
these  time  histories  at  a time  the  spin  was  considered  to  be  in  a steady,  flat  mode.  This 
time  was  selected  at  85  sec  into  the  spin. 

Flight  motions  calculated  using  measured  static  and  forced-oscillation  aerodynamic 
data  (Aero  Model  1)  are  shown  in  comparison  with  the  flight  test  motion  in  Fig.  3.  The 
technique  used  in  extracting  the  flight  test  components  is  presented  in  Appendix  C.  It 
should  be  noted  here  that  the  symbols  presented  for  the  calculated  characteristics  are 
cosmetic,  whereas  the  solid  line  represents  the  actual  calculated  points.  As  can  be  seen  in 
this  figure,  the  calculated  motions  generated  require  approximately  40  sec  to  stabilize. 
The  motions  beyond  40  sec  can  be  considered  representative  of  the  flight  test  spin 
characteristics.  The  perturbations  in  the  calculated  motions  during  the  first  40  sec  can  be 
attributed  to  the  lack  of  balance  between  the  aerodynamic,  inertia,  and  thrust 
contributions  in  the  moment  equations  of  motion.  The  unsteadiness  can  be  seen  in  Fig. 
3,  where  the  aerodynamic  contributions  Paero>  <1a ero>  and  t'aero  to  the  total 
moment  equations  of  motion  are  presented  along  with  flight  test  data. 

Flight  motions  calculated  using  measured  static  and  forced-oscillation  aerodynamic 
data,  in  which  the  forced-oscillation  derivatives  were  restricted  to  the  oscillatory 
components  of  the  total  angular  rates  (Aero  Model  2),  are  shown  in  Fig.  4.  The  flight 
motions  are  also  shown  for  comparison.  The  overall  spin  characteristics  are  changed  by 
Aero  Model  2 and  result  in  a gradual  reduction  in  angle  of  attack  and  total  angular  rate 
(J2)  with  time.  The  motion  is  suppressed  to  a nonoscillatory  state  which  is  not 
characteristic  of  the  flight  test  spin.  The  suppression  of  Paero>  Qaero,  and  r a e r o to 
a nonoscillatory'  state  strongly  indicates  the  importance  of  including  the  influence  of 
steady  rotation  on  the  aerodynamics. 

Flight  motions  were  calculated  by  incorporating  the  rotation-balance  aerodynamic 
data  (Aero  Model  3).  A comparison  of  the  calculated  and  flight  test  motions  is  shown  in 
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Fig.  5.  Although  the  attitude  angles  and  rotation  rates  of  the  calculated  motions  are  not 
exact  overlays  of  the  flight  test  motions,  these  results  are  considered  quite  representative 
since  completely  unaltered  aerodynamic  data  were  used.  A brief  examination  of  the  time 
history  reveals  that  the  aerodynamic  modeling  as  shown  by  the  calculated  terms  Paf.ro. 
4a fro.  and  ^AERO  more  closely  duplicates  the  flight  characteristics  than  either  Aero 
Model  1 or  2.  Some  discrepancy  in  the  longitudinal  aerodynamic  model  is  indicated  by 
the  lack  of  correlation  between  the  calculated  and  flight  test  values  of  cjaero- 

The  general  conclusions  previously  drawn  indicate  that  for  the  specific  flight  test 
motions  simulated,  use  of  either  Aero  Model  1 or  3 would  yield  representative  spin 
motions.  This  general  conclusion  resulted  in  an  attempt  to  answer  the  question  whether 
representative  spin  motions  could  be  obtained  using  Aero  Models  1 and  3 for  other  initial 
conditions.  The  thought  here  was  to  introduce  a change  in  the  initial  conditions  that 
would  yield  a known  result.  This  led  to  the  application  of  left  pro-spin  controls  (LPSC) 
at  the  start  of  the  calculated  motion.  In  general,  application  of  pro-spin  controls  should 
result  in  an  increase  in  the  spin  rate  and  angle  of  attack  with  some  reduction  in  the 
oscillatory  characteristics. 

Figure  6 presents  a comparison  of  the  calculated  motions  using  Aero  Models  1 and  3 
for  application  of  LPSC.  For  Aero  Model  1.  the  application  of  LPSC  resulted  in  driving 
the  simulated  airplane  out  of  the  spin.  This  result  can  be  explained  by  the  control  power 
characteristics  in  the  50-  to  65-deg  angle-of-attack  range;  both  rudder  and  aileron  produce 
a retarding  moment  in  this  angle-of-attack  range.  The  initial  perturbations  resulting  from 
application  of  LPSC  drives  the  calculated  motion  into  this  angle-of-attack  region. 
Application  of  LPSC  to  Aero  Model  3 increased  the  spin  rate.  This  influence  resulted  in  a 
gradual  suppression  of  the  oscillatory  characteristics  along  with  an  increase  in  angle  of 
attack  and  spin  rate  (£2).  The  resultant  motions  for  Aero  Model  3 were  as  expected  for 
application  of  pro-spin  controls,  whereas  for  Aero  Model  1 they  were  not.  Application  of 
pro-spin  controls  to  Aero  Model  3 improved  the  correlation  with  the  flight  test  spin.  This 
improvement  in  the  flight  correlation  resulting  from  application  of  LPSC  indicates 
possible  discrepancies  in  the  static  aerodynamic  model. 

An  attempt  was  made  to  simulate  the  steady,  flat  spin  using  Aero  Model  3 without 
the  forced-oscillation  derivatives.  A comparison  of  the  calculated  motions  with  and 
without  forced-oscillation  derivatives  is  shown  in  Fig.  7.  The  calculated  motions  without 
forced-oscillation  derivatives  were  terminated  at  8 sec  because  of  the  oscillatory 
divergence  that  resulted  in  exceeding  the  aerodynamic  data  range  in  sideslip.  This  result 
indicates  the  necessity  of  including  some  type  of  damping  data  to  handle  oscillations  in  p 
arid  q.  Oscillatory  damping  data  for  both  pitch  and  roll  could  not  be  obtained  by  the 
technique  used  in  acquiring  the  steady-state  rotation-balance  data. 
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Further  investigation  was  made  to  determine  which  of  the  incremental  moment 
contributions  (ACgrb , ACni  ib , or  AC„rb)  caused  by  steady-state  rotation  were  most  important 
in  maintaining  the  steady  spin.  Figures  8,  9,  and  10  present  the  results  obtained  for 
steady  spin  attempts  using  Aero  Model  3 when  these  moment  contributions  are 
implemented  separately.  The  calculated  motion,  which  includes  all  force  and  moment 
contributions  caused  by  steady  rotation,  shown  earlier  in  Fig.  5,  is  presented  in  Figs.  8, 
9,  and  10  for  comparison.  An  examination  of  Figs.  8 and  10  reveals  that  the  results  for 
individual  implementation  of  the  rolling  and  yawing  moment  contributions  (ACglb  and 
ACnrb)  caused  by  steady  rotation  are  quite  representative  of  those  including  all 
rotation-balance  data.  In  general,  the  oscillatory  characteristics  are  slightly  more 
suppressed  for  the  inclusion  of  ACnrb  than  for  ACgrb  alone.  The  calculated  motions 
using  Aero  Model  3 for  ACmib  alone  result  in  a gradual  reduction  in  angle  of  attack  and 
total  angular  rate  (£2).  The  motion  is  suppressed  to  a non  oscillatory  state  which  is  not 
representative  of  the  calculated  motions,  including  all  rotation-balance  data  as  well  as  the 
flight  test  spin.  A comparison  of  the  calculated  motions  presented  in  Figs.  4 and  9 
confirms  that  the  pitching-moment  contribution  caused  by  steady  rotation,  ACmrb,  alone 
has  little  influence  on  the  calculated  motions.  For  this  configuration,  the  most  significant 
factor  that  must  be  included  in  the  aerodynamic  model  for  steady  spin  calculation  is  the 
nonlinear  variation  of  the  rolling-  and  yawing-moment  contributions  caused  by  rotation 
rate. 


Once  it  was  established  that  the  rotation-balance  data  system  (Aero  Model  3)  was 
best  for  calculating  the  steady,  flat  spin  for  the  test  configuration,  an  attempt  was  made 
to  determine  possible  sources  of  error  preventing  better  calculated/flight  test  correlation. 
Previous  conclusions  indicated  that  the  most  obvious  discrepancy  was  in  the  longitudinal 
aerodynamic  model.  This  is  better  pointed  out  by  the  lack  of  correlation  between  the 
calculated  and  flight  test  values  of  qAERO  shown  in  Fig.  5.  The  incremental  difference 
between  mean  values  of  flight  test  and  calculated  4aero  corresponds  to  an  increment  in 
pitching  moment  between  0.12  and  0.17. 

Results  from  the  application  of  LPSC  in  the  calculated  motions  indicated  possible 
discrepancies  in  the  static  aerodynamic  data.  Therefore,  additional  static  longitudinal 
wind  tunnel  data  (normal-force  and  pitching-moment  data)  were  obtained  from  a 
1/15-scale  model  F-4D  test  at  NASA  Ames  (Ref.  13).  The  NASA  Langley  data  used  in 
the  previously  calculated  motions  are  presented  along  with  the  NASA  Ames  data  in  Figs. 
11  and  12  for  comparison.  The  NASA  Ames  test  was  conducted  at  a Reynolds  number 
of  4.0  x 106  per  foot,  whereas  the  Langley  test  was  conducted  at  a Reynolds  number  of 
1.35  x 106  per  foot.  Assuming  that  the  differences  between  the  Ames  and  Langley  data 
are  not  caused  by  model  configuration,  the  difference  in  the  normal-force  and 
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pitching-moment  data  at  angles  of  attack  above  approximately  20  deg  (stall)  is  most 
probably  caused  by  Reynolds  number  and/or  sting/support  interference,  as  indicated  in 
Ref.  14.  A comparison  of  the  incremental  differences  between  the  Ames  and  Langley 
data  in  the  stall  angle-of-attack  region  gives  credence  to  the  possibility  of  Reynolds 
number  effects.  In  the  higher  Reynolds  number  test,  an  increase  in  lift  occurs  with  a 
resulting  reduction  in  pitching  moment.  In  the  spin  angle-of-attack  region  (a  = 55  to  75 
deg),  the  increments  in  normal  force  and  pitching  moment  vary  between  -0.10  and  -0.15 
for  ACZ  and  between  0.125  and  0.175  for  ACm.  Somewhat  of  a surprise  is  the  fact  that 
the  incremental  difference  betwen  the  mean  values  of  flight  test  and  calculated 
aerodynamic  pitching  moment  (Fig.  5)  corresponds  to  the  increments  between  the  Ames 
and  Langley  data  (Fig.  12).  It  should  be  noted  that  the  average  Reynolds  number  for  the 
flight  test  spin  was  14.0  x 106. 

Calculations  were  made  to  determine  the  effects  on  the  steady  spin  motion  of  these 
increments  in  static  normal  force  and  pitching  moment.  Spin  motions  were  calculated 
using  Aero  Model  3 and  the  minimum  and  maximum  values  of  the  increments  in  Cz  and 
Cm  between  the  Langley  and  Ames  data.  Figure  13  presents  a comparison  of  the  spin 
calculations  for  the  minimum  and  maximum  values  of  the  increments  in  Cz  and  Cm 
between  the  Langley  and  Ames  data.  The  overall  results  for  the  calculated  spins  for  each 
set  of  increments  are  not.  significantly  different.  Figure  14  presents  a comparison  of 
calculated  spins  for  Aero  Model  3 and  Aero  Model  3 with  increments  in  ACZ  of  -0.15 
and  in  ACm  of  0.175.  The  flight  test  motion  is  also  shown  in  Fig.  14  as  a means  of 
determining  whether  any  improvements  in  correlation  were  made.  The  results  obtained 
for  increments  of  ACZ  = -0.15  and  ACm  = 0.175  significantly  improve  the  flight  test 
correlation  and  strengthen  the  basis  for  requiring  model  data  corrected  for  such  effects  as 
Reynolds  number,  sting  and  wall  interference,  and  model  configuration  differences. 

The  above  results  indicate  that  a good  static  data  matrix  is  necessary  for  better 
correlation  with  flight  test  data,  but  does  not  eliminate  the  predominant  oscillatory 
characteristic  in  pitch.  An  attempt  at  determining  the  source  of  error  led  to  a review  of 
the  flight  test  aerodynamic  data.  In  particular,  the  contributions  to  the  flight  pitch 
angular  acceleration  due  to  thrust,  aerodynamics,  and  inertial  coupling  determined  from 
the  time  history  are  shown  in  Fig.  15.  The  oscillatory  characteristic  of  flight  test 
Qtotal  is  dramatically  influenced  by  the  inertial  coupling  component  and  not  by  the 
aerodynamic  component.  This  indicates  that  the  oscillatory  characteristic  of  the 
calculated  qAERO  as  shown  in  Fig.  14  could  be  caused  by  a discrepancy  in  the 
forced-oscillation  pitch-damping  wind  tunnel  data.  In  light  of  this  possibility,  a review  of 
the  forced-oscillation  pitch  data  (solid-line  fairing  in  Fig.  16)  used  in  the  calculations 
revealed  a large  reduction  in  (Cmq  + Cm-a)  in  the  spin  angle-of-attack  region  (a  = 55  to 
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80  deg).  However,  additional  forced-oscillation  data  (long/short  dash  fairing)  obtained  on 
a 1/11-scale  F-4B  model  at  a Reynolds  number  of  0.45  x 102 * * * 6  per  foot  did  not  exhibit 
the  same  trend  in  the  spin  angle-of-attack  region;  a comparison  of  the  1/1 1-scale  data 
with  that  used  for  the  calculated  motion  is  shown  in  Fig.  16.  The  agreement  in  the  data 
up  to  approximately  stall  angle  of  attack  (a  « 23  deg)  indicates  that  differences  in  the 
stall/spin  angle-of-attack  region  are  not  caused  by  model  configuration  but  possibly  are 
caused  by  Reynolds  number,  sting,  and  wall  interference  effects.  The  incremental 
differences  in  the  stalled  angle-of-attack  region  are  characteristics  of  Reynolds  number 
effect.  Since  the  0.1 3-scale  data  were  obtained  at  a higher  Reynolds  number,  they  were 
retained  and  modified  as  indicated  by  the  dashed-line  fairing  in  the  55-  to  80-deg 
angle-of-attack  region. 

Figure  17  presents  a comparison  of  the  time  history  spin  calculations  using  Aero 
Model  3 with  and  without  modifications  in  the  static  and  forced-oscillation  pitch  data 
[ACZ  = -0.10,  ACm  = 0.125,  (Cmq  + Cmi)  fairing]  and  the  flight  test  time  history.  The 
resultant  motions  for  the  modified  data  more  closely  model  the  longitudinal  flight 
characteristics.  This  is  shown  in  Fig.  17  by  the  much  improved  relationship  between  the 
calculated  and  flight  test  values  of  a,  q,  and  qAERO-  It  should  be  pointed  out  that  with 
these  modifications,  the  calculated  motion  exhibits  a slight  overdamping  in  the 
lateral/directional  characteristics.  The  indications  are  that  possible  modifications  in  the 
lateral/directional  data  matrix  would  be  necessary.  In  general,  the  overall  calculated/flight 
test  correlation  is  much  improved  by  the  data  matrix  modifications. 

4.0  CONCLUDING  REMARKS 

The  results  of  this  analytical  investigation  of  steady  spin  characteristics  of  the  F-4 
series  airplane  using  conventional  and  rotation-balance  aerodynamic  data  are  summarized 
as  follows. 

1.  For  the  configuration  investigated,  the  use  of  conventional  static  and 
forced-oscillation  aerodynamic  data  as  combined  in  Aero  Models  1 and  2 
would  not  result  in  maintaining  the  steady,  flat  spin  motions. 

2.  The  inclusion  of  rotation-balance  data  in  the  aerodynamic  data  matrix 

(Aero  Model  3)  resulted  in  calculated  motions  that  are  considered 

representative  of  the  flight  test  steady  spin.  The  inclusion  of 
forced-oscillation  data  to  account  for  the  effects  of  oscillations  in  roll  and 
pitch  superimposed  on  the  steady  spin  motion  were  found  to  be  necessary 

in  order  to  maintain  the  steady  spin.  For  this  configuration,  the  most 

significant  factor  that  must  be  included  in  aerodynamic  modeling  for 
steady  spin  simulation  is  the  nonlinear  variation  of  the  rolling  and  yawing 
moment  contributions  caused  by  steady  rotation. 
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3.  The  major  deficiency  in  utilizing  rotation-balance  data  in  this  type  of 
simulation  is  the  lack  of  a complete  set  of  measured  data.  The 
rotation-balance  data  at  sideslip  angles  up  to  15  deg  were  quite  adequate 
for  the  flight  test  spin  investigated;  however,  measured  rotation-balance 
data  should  be  obtained  at  larger  angles  of  sideslip  for  application  in  the 
flight  regime  between  initial  stall  and  steady  spins.  Rotation-balance  data 
should  be  obtained  down  to  at  least  stall  angle  of  attack  for  attempted 
simulation  of  stall/post-stall  gyration  maneuvers. 
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SPIN  AXIS 


Figure  2.  Rotation-balance  apparatus. 
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Figure  3.  Comparison  of  calculated  motions  using  Aero  Model  1 
and  airplane  flight  motions. 


19 


AEDC-T  R -7  7-1 2 6 


Ti 

i/i 

S' 

■o 


O 

Q£ 


■ 

s 

ss 

YVi 

<Vk 

AA 

A/1 

50 


S 

V) 

S’ 


o 

Of 


-50 


,dkL 

mu 

B— 

H 

\0'^qjfcrLv 

WHr«a  Lviaa 

rafe^SI 

10  15  20  25  30  35  40 

Time,  sec 


50 


Figure  3.  Concluded. 
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Figure  4.  Comparison  of  calculated  motions  using  Aero  Model  2 
and  airplane  flight  motions. 
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Figure  5.  Comparison  of  calculated  motions  using  Aero  Model  3 
and  airplane  flight  motions. 
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Figure  6.  Comparison  of  calculated  motions  for  Aero  Models  1 
and  3,  left  pro-spin  controls. 
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o Aero  Model  3 - Without  Damping 
4 Aero  Model  3 - With  Dampl  ng 
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Figure  7.  Comparison  of  calculated  motions  using  Aero  Model  3 
with  and  without  forced -oscillation  damping  data. 
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Figure  9.  Comparison  of  calculated  motions  using  Aero  Model  3 
system  with  ACm  rb  only  and  rotation-balance  data. 
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Figure  10.  Comparison  of  calculated  motion  using  Aero  Model  3 
system  with  AC„rb  only  and  rotation^balance  data. 
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Figure  10.  Concluded. 
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Figure  11.  Comparison  of  NASA  Ames  and  NASA  Langley  static  normal-force 
coefficient,  controls  neutral  position. 
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Figure  12.  Comparison  of  NASA  Ames  and  NASA  Langley  static  pitching-moment 
coefficient,  controls  neutral  position. 
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Figure  13.  Comparison  of  calculated  motions  using  Aero  Model  3 for 
longitudinal  coefficient  increments. 


37 


33S/PEJ  '0  tap 


• deg ,sec  q,  deg/sec  p,  deg /sec 


AEDC-TR-77-1  26 


0 5 10  15  20  25  30  35  40  45  50 

Time,  sec 

Figure  14.  Comparison  of  calculated  motions  using  Aero  Model  3 with 
and  without  longitudinal  increments  and  flight  test. 
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o Aero  Model  3 - ACZ  ■ -0. 15,  ACm  ■ 0. 175 
a Aero  Model  3 
Flight  Test 


Figure  14.  Concluded. 
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Figure  15.  Components  of  flight  test  pitch  angular  acceleration. 
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Figure  16.  Variation  of  aerodynamic  damping  in  pitch  with  angle  of  attack. 
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Figure  17.  Comparison  of  calculated  motions  using  Aero  Model  3 
with  and  without  modifications  in  the  static  and 
pitch-damping  data  and  flight  test. 
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Figure  17.  Concluded. 
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Table  1.  Mass,  Inertia,  and  Dimensional  Characteristics  of  the  F-4  aircraft 


Loading  - AIM-7' s on  Armament  Stations  3 and  7 

Weight,  lb 

eg,  percent  MAC  

Moments  of  Inertia,  slug-ft2 


Wing: 

Span,  ft  

Area,  ft^  

Mean  Aerodynamic  Chord,  ft 


Inertia  Coupling  terms: 


*Y  " rz 


rz  ~ *X 


rx  “ JY 


-0.836 


0.977 


= -0.772 


*XZ 

h 

^XZ 

*xz 


0.241 


0.033 


0.030 


36,950 

33.7 

22,600 

163,100 

182,000 

5,450 


38.41 

538.34 

16.04 
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Table  2.  Initial  Conditions  from  Flight  202,  Run  12  of  Ref.  12 


Conditions  at  85 

= 74  deg 
0=1  deg 
9 = -28.0  deg 
<t>  =8.0  deg 
^ = 0 deg 
q = -20.0  deg/sec 
p = -30.0  deg/sec 
r = -82.0  deg/sec 
q = 15.0  deg/sec2 


>c  in  the  Spin: 

P 

r = 

h 

q 

Moo  = 

d 

T 

€ = 5.25  deg 


-62.0  deg/sec2 
—11.0  deg/sec2 
12,500  ft 
35  lb/ft2 
0.18 

-0.3866  ft 
62  lb  (at  Idle  rpm) 
42  slug-ft2 
518.0  rad/sec 


Controls  Time  Function: 


Time,  sec 

0.0 

n 

6h,  deg 

6r.  deg 

6a,  deg 

3 0 

■ ou  • u 

0 0 

3 .0 

0.0  

► 

0.0 

0.0 
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A.1  SIX-DEGREE-OF-FREEDOM,  RIGID-BODY  EQUATION  OF  MOTIONS 


The  dynamic  equations  required  to  specify  the  translational  and  rotational  motions 
of  a rigid  body  moving  through  space  are  described  in  this  appendix.  The 
six-degree-of-freedom  nonlinear  differential  equations  representing  linear  and  angular 
accelerations  on  a moving  body  axis  system  (Fig.  A-l)  having  its  origin  at  the  aircraft 
center  of  mass  are  presented  below. 


Forces: 


. „ f r 

u = rv  - qw  - e sin  6 + _*  - 

m rr. 


v = pw  - ru  + g cos  d sin  <£  + 


U 

m 


, F T 

w = qu  — pv  + g cos  u cos  <p  + — — + —5. 


Moments: 


p = 1.1 [L  qr  + 2 (f  + pq)  + HjL 

Ix  Ix  Ix 


q = 'z~'x  pr  + _x_z^  _ p^>  + 21  + !ILXX  + 1,1  VG 


T = lx lx  p q + *XZ  (p  _ qr)  + !1z  t H^ZC. 

1 2 lz  tz  Iz 


The  external  forces  and  moments  (Fx,  FV)  Fz,  Mx>  My,  and  Mz)  in  the  equations  are 
comprised  of  aerodynamic  coefficients  representative  of  the  aircraft.  The  external  force 
and  moment  contributions  due  to  engine  thrust  (including  gyroscopic  effect)  are 
represented  by  TX)  Tz>  Myt>  Myg.  ar>d  MZg  and  are  developed  in  Appendix  A.2. 
Development  of  the  aerodynamic  math  model  used  is  presented  in  Appendix  A.3.  The 
method  used  in  implementing  the  rotational-balance  data  is  outlined  in  Appendix  B. 

Auxiliary  equations  used  in  this  analysis  are  given  as  follows: 
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0 = 


p - tan  6 (r  cos  0 + q sin  0) 


* = 


r cos  0 + q sin  0 
cos  Q 


O _ up  + vq  + wr 

ss  " N 

An  explanation  of  the  steady-state  rotational  rate  vector  (^ss)  is  given  in  Appendix  B.2. 

The  equations  of  motion  are  numerically  intergrated  to  provide  time  histories  of  the 
aircraft  spinning  motion. 

A.2  EQUATIONS  DEFINING  THE  FORCE  AND  MOMENT  CONTRIBUTIONS  DUE  TO 
ENGINE  THRUST 


deg 

Line 
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Forces 


Moments 


Tx  = T cos  f 
l'Y  = 0 

T y = -rr  sin  ( 


MXT  = 0 
My  1-  = T(d) 

Myg  = -*xh^Er 
V,ZG  - 


where 


<f>  Y = Engine  rotor  angular  velocity,  rad/sec 

^XE  = Engine  rotor  moment  of  inertia,  slug/ft2 
d = Distance  between  thrust  line  and  eg 

Assumptions 

1.  Engine  thrust  line  parallel  to  X-Z  plane. 


,2.  . Gyroscopic  moments  assume  engine  axis  parallel  to  X axis. 

3.  No  variation  in  Y and  Z eg  location. 

A.3  EQUATIONS  DEFINING  THE  TOTAL  AERODYNAMIC  DATA  ALONG  AND 
ABOUT  EACH  BODY  AXIS 

LONGITUDINAL  AXIS  PLANE 


Fx  . qS  [c>,0,  8„)  + AC,rt(.,  0,<5H.  %i)  + (c%(a))il] 
Fz  - qS  [c.  (a.  0,  S„>  + AC  ^(a,  0,  *„ , + (c,w)  !£] 

My  . qSC  [C.U,  0,  S„)  + AC,,,^  («,  JS,*,,.  Sat)  + (c,%<«>) 
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LAT  ERA I^DI  RECTION AL  AXIS  PLANE 


F 


Y 


+ 


2V  / 


+ 


Mz  = qSb  |Cn(a,  /3)  + ACn(a,  /S , 8a)  * ACn(a,  8) 


AC^H’¥MVa)]  t£ + tV-M 


+ c . + c. 

n/3  2 V 


q C 

o 

‘q  2 V 


My  = qSb  JC.  (a,  /9)  + AC.  (a,  jS  , SB)  + AC  Co,  /S  , 8) 


r b 
_o 

2 V 


Cn  b + C ‘1oG  1 

X Tv”  i 


where 


ACgCo,  j8,  Sa)  = ACg(o,  /3,  30) 
right  roll  command. 


surface  command  determines  whether  left  or 


ACg(a,  /S,  8)  = ACg  (a,  0 , 1 30)  (jr) 


where  the  "A"  equations  are  the  same  as  for  Cg  above  except  Cg  is  replaced  by  Cy  or  C„. 
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A.4  EQUATIONS  FOR  TRANSFERRING  AERODYNAMIC  DATA  INPUTS  TO  THE 
PROPER  HORIZONTAL  CENTER  OF  GRAVITY  POSITION 


where 


Mx  = Mx 

eg 

My  = My  - Fz  (RI)ELCG) 

eg 

Mz  = Mz  + Fy  (RDELCG) 


RDEjLCG  = — cSnn’.r^  ^ ^ 


Figure  A-l.  Body  axis  system. 
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APPENDIX  B 

IMPLEMENTATION  METHODS  FOR  ROTATION-BALANCE  DATA 

B.1  ADDITION  OF  INCREMENTAL  ROTATION-BALANCE  DATA  TO  AERODYNAMIC 
FORCE  AND  MOMENT  EQUATION 


AC. 


rb 


(■•  tt)  ■ 4C*„,  (■•  4fr  [4Sb  (-  V) 


- 4C-,b  (■•  *• 


where 


4C»,b(“^'-21’^)- 


incremental  rotation— balance  data  for  longitudinal  control 
surface  neutral 

incremental  rotation— balance  data  for  longitudinal  control 
surface  at  -21  deg 


For  other  coefficients,  replace  Cm  with  Cx,  Cz,  CY>  C„,  or  Q. 
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B.2  EQUATIONS  DEFINING  THE  STEADY-STATE  ROTATIONAL  RATE  VECTOR 

The  rotation-balance  data  incorporated  in  the  aerodynamic  data  model  of  Appendix 
A.3  were  obtained  with  a model/balance  rotating  in  the  wing  tunnel  with  its  rotational 
vector 'aligned  with  the  velocity  vector  of  the  wind  tunnel.  Because  the  rotational  and 
velocity  vectors  are  normally  not  aligned  in  spinning  flight,  the  direct  implementation  of 
the  rotation  balance  data  as  a function  of  the  total  rotational  rate  vector  (£2)  is  in  error. 
The  method  used  for  implementing  the  data  resolves  a component  of  the  total  rotational 
vector  along  the  aircraft  velocity  vector.  This  component  rotational  vector  (£2SS)  may 
now  be  used  in  determining  the  incremental  coefficients  to  be  included  in  the 
aerodynamic  data_jnatrices.  The  determination  of  the  magnitude  of  the  steady-state 
rotational  vector  (£2SS)  and  its  inclination  (6)  to  the  aircraft  total  rotational  vector  (£2)  is 
outlined  below  and  is  the  same  as  that  used  in  Ref.  15. 

X Body  Axis 


Y Body  Axis 

J « 


v = ul  + v]  + wlc  = Total  velocity  vector 

£l=pi  + qj  + rE  = Total  rotational  vector 

nss  = pssT  + qsy/  + rssE  = Steady-state  rotational^vector 

along  velocity  vector  V. 

Ift  I = Ifil  COS  5,  COS  8 = i^ssl 

! &s  1 ' |Q| 

£2  • V = 1 0 1 | V | cos  8 

up  - vq  + wr  = jilj  jV'l  cos  5 
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Now  substituting. 


But 


If)  | __  up  4-  vq  + wr 

SS|  || 

|V| 


JL  = cos  B cos  a 

V 

$ = sin  P 

= cos  B sin  a 
* 


SS 


u 


X Body 
Axis 


Y Body 
Axis 


again  substituting, 


«ss  = p cos  y3  cos  a + q sin  j8  + r cos  £ sin  a 


_^To  separate  p,  q,  and  r into  steady-state  (parallel  to  V)  and  oscillatory  (perpendicular 
to  V)  components,  the  components  of  £2SS  along  the  X,  Y,  and  Z axes  are  first  determined. 

i = i \ 

= Qss  c°s  p c°s  a i + flaa  sin 

- fiss  cos  yS  sin  a E 

Components  of  fiss  in  thetT  and  k directions  are  the  steady-state  components  of  p,  q,  and 
r,  respectively. 

Therefore 

Pss  = Qss  cos  P cos  « 

<Us  = ftss  sin  P 

rss  = Qss  cos  P sin  a 
The  oscillatory  components  are: 

Po  = P - Pss 

= q - %s 


r0  = r - rss 
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The  oscillatory  components  are  used  in  conjunction  with  the  dynamic  oscillatory 
derivatives  in  the  aerodynamic  data  matrices  of  Appendix  A.3.  Utilizing  the  rotational  and 
oscillatory  derivatives  in  the  manner  outlined  accounts  for  oscillations  that  are 

■ ■ .it . ' 

superimposed  on  an  aircraft  steady-state  spinning  motion. 

In  the  subject  investigation,  the  rotation-balance  data  were  applied  when  the 
inclination  angle  (8)  between  the  rotational  vector  (£2)  and  velocity  vector  (V)  was  less  than 
40  deg.  Above  40  deg,  the  rotation-balance  increments  to  the  aircraft  aerodynamic  matrix 
were  set  equal  to  zero  with  an  exponential  decay  function  (see  Appendix  B.3).  Under  these 
conditions,  the  oscillatory  rates  (p0,  q0,  and  r0)  become  equal  to  the  total  rotational  rates 
(p.  q,  and  r).  Likewise,  when  the  total  rotational  vector  (i2)  is  aligned  with  the  velocity 
vector  (V),  the  condition  where  5 = 0,  the  aircraft  is  in  a flat  spin  and  the  oscillatory  rate 
components  are  zero. 

B.3  ELIMINATION  OF  DISCONTINUITY  EFFECTS  WHEN  IMPLEMENTING 
ROTATION-BALANCE  DATA 


. A problem  that  arises  in  using  the  inclination  angle  criteria  is  the  need  to  minimize  the 
discontinuity  which  occurs  when  the  criteria  for  use  of  rotary  data  are  satisfied  when  large 
rotation  rates  exist.  At  such  points,  significant  magnitudes  of  incremental  rotation-balance 
aerodynamics  can  be  suddenly  added  into  the  equation  as  the  rotation  vector  swings  inside 
of  the  limiting  angle.  The  addition  of  this  data  is  instantaneous  and  is  therefore  not  realistic. 
To  minimize  this  discontinuity,  an  exponential  decay  was  used  to  smoothly  bring  the 
incremental  rotation-balance  data  into  the  equations  of  motion.  By  using  a small  time 
constant,  the  exponential  function  eliminates  any  pulse  effects  yet  does  not  reduce  the  total 
effect  of  the  rotation-balance  data.  The  technique  is  outlined  below. 


AC 


'rb 


AC 


- AC. 


_ — Ai/j 


rb 


sm  j 


where 


AC. 


AC. 


rb  j 


AC. 


s m j—  1 


fj  = time  iteration  at  which  enter  or  exit 


I 


rotation— balance  data  At 


li) 


iteration  at  previous  entry  or  exit 


t = ftashout  time  constant 


AC.. 


Mi 


(is  defined  in  Appendix  B.l) 


For  other  moment  coefficients,  replace  "Cm  " with  "Cg"  or  "Cn."  The  above  technique  was 
developed  and  presented  in  Ref.  1 5. 
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APPENDIX  C 

DETERMINATION  OF  AERODYNAMIC  MOMENTS  FROM  FLIGHT  TEST 


In  order  to  determine  angular  accelerations  produced  only  by  aerodynamics, 

M v 

Pafro  = f- 

X 


My 

qAFRO  _ T- 
*Y 


M ? 

rAFRO  " 

1 Z 


the  inertial  and  engine  effects  were  subtracted  from  the  calculated  body  axis  accelerations 
(slope  of  measured  body  axis  angular  rates). 


Paf.ro  = P “ v - d + pq) 


'»*«o  ' ^ ^ ^ - T7  K * "t.) 


fAF.RO  ‘ ! - H - ¥ (f  - * - f (MZr) 

z 'z  *z  v G f 


The  pitch-,  roll-,  and  yawing-moment  coefficients  attributable  to  all  aerodynamic  terms  are 


(C  ) = qAEHO  |y 

,n  AERO  q SC 

(C0)  = Paero  *X 

* AERO  q"  Sb 


(C  ) = fAERO  !Z 

nAERO  q Sli 
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NOMENCLATURE 

vr 

b 

Wingspan,  ft 

C 

Wing  mean  aerodynamic  chord,  ft 

CGNEW 

New  center-of-gravity  location,  percent  chord 

CGREF 

Reference  center-of-gravity  location,  percent  chord 

Q 

Rolling-moment  coefficient,  rolling  moment/qSb  about  airplane  eg 

C‘p 

Derivative  of  rolling-moment  coefficient  with 
9Q2/d(pb/2V),  per  radian 

respect  to  roll 

rate. 

Derivative  of  rolling-moment  coefficient  with 
3Cg/3(qC/2V),  per  radian 

respect  to  pitch 

rate, 

cer 

Derivative  of  rolling-moment  coefficient  with 
3Cg/3(rb/2V),  per  radian 

respect  to  yaw 

rate, 

Derivative  of  rolling-moment  coefficient  with  respect  to  0,  3C£/3(j3b/2V),  per 

radian 


Cm  Pitching-moment  coefficient,  pitching  moment/qSC  about  airplane  eg 

Cmp  Derivative  of  pitching-moment  coefficient  with  respect  to  roll  rate, 

9Cm/3(pb/2V),  per  radian 

Cm'q  Derivative  of  pitching-moment  coefficient  with  respect  to  pitch  rate, 

/3(qC/2V),  per  radian 

Cm  r Derivative  of  pitching-moment  coefficient  with  respect  to  yaw  rate, 

3Cm/3(rb/2V),  per  radian 

Cm  -a  Derivative  of  pitching-moment  coefficient  with  respect  to  a,  3Cm/3(aC/2V), 

per  radian 

Cn  Yawing-moment  coefficient,  yawing  moment/qSb  about  airplane  eg 

Cnp  Derivative  of  yawing-moment  coefficient  with  respect  to  roll  rate, 

3Cn/3(pb/2V),  per  radian 

C„q  Derivative  of  yawing-moment  coefficient  with  respect  to  pitch  rate, 

3Cn/3(qC/2V),  per  radian 
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Cnt  Derivative  of  yawing-moment  coefficient  with  respect  to  yaw  rate, 

3Cn/3(rb/2V),  per  radian 

C„£  Derivative  of  yawing-moment  coefficient  with  respect  to  0,  3C„/3(0b/2V), 

per  radian 

Cx  Longitudinal-force  coefficient,  longitudinal  force/qS 

CXq  Derivative  of  longitudinal-force  coefficient  with  respect  to  pitch  rate, 

3Cx/3(qC/2V),  per  radian 

CX(j  Derivative  of  longitudinal-force  coefficient  with  respect  to  a,  3Cx/3(aC/2V), 

per  radian 

Cy  Side-force  coefficient,  side  force/qS 

Cy  p Derivative  of  side-force  coefficient  with  respect  to  roll  rate,  3Cy /3(pb/2V), 

per  radian 

CVl  Derivative  of  side-force  coefficient  with  respect  to  yaw  rate,  3Cy /3(rb/2V), 

per  radian 

' 

Cyp  Derivative  of  side-force  coefficient  with  respect  to  0,  3Cy /3(0b/2V),  per 

. radian 

Cz  Normal-force  coefficient,  normal  force/qS 

CZq  Derivative  of  normal-force  coefficient  with  respect  to  pitch  rate, 

3Cz/3(qC/2V),  per  radian 

CZ(j  Derivative  of  normal-force  coefficient  with  respect  to  a,  3Cz/3(aC/2V),  per 

radian 

eg  Center-of-gravity  location,  percent  chord 

d Distance  from  thrust  line  to  center  of  gravity 

Fx  Force  acting  along  X-body  axis,  lb 

Fy  Force  acting  along  Y-body  axis,  lb 

Fz  Force  acting  along  Z-body  axis,  lb 

g Acceleration  of  gravity,  ft/sec2 
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h.. 

!x  Jy  ,h 

Ixz 

IXE 

ij.k  .•>  . 
k 

LPSC 

M 

oo 

Mx  ■ , 

My 

. / 

> ' 1 . 

Myg 

My  t 

Mz 

Mzg 

f ■ 

m 

P,Q,r.  t 
Paero 
Qaero 
taero 
Po  :Qo  ,r0 

Pss>qss>rss 

q 


Altitude,,  ft 

Moments  of  inertia  about  X-,  Y-,  and  Z-body  axes,  respectively,  slug-ft2 
Product  of  inertia,  slug-ft2 

Moment  of  inertia  about  X-body  axis  due  to  engine  rotation,  slug-ft2 
Unit  vectors  along  X-,  Y-,  and  Z-body  axes,  respectively 
Frequency  parameter,  fib/2V  or  flC/2V 
Left  pro-spin  controls 
Mach  number 

Moment  acting  about  X-body  axis,  ft-lb 
Moment  acting  about  Y-body  axis,  ft-lb 

Moment  acting  about  Y-body  axis  caused  by  engine  rotation,  ft-lb 
Moment  acting  about  Y-body  axis  caused  by  engine  thrust,  ft-lb 
Moment  acting  about  Z-body  axis,  ft-lb 

Moment  acting  about  Z-body  axis  caused  by  engine  thrust,  ft-lb 
Mass,  slugs 

Components  of  ^ about  X-,  Y-,  and  Z-body  axes,  respectively,  rad/sec 
Aerodynamic  contribution  to  p,  deg/sec2 
Aerodynamic  contribution  to  q,  deg/sec2 
Aerodynamic  contribution  to  r,  deg/sec2 

Components  of  S2  about  X-,  Y-.  and  Z-body  axes  caused  by  oscillations 
superimposed  on  steady  rotation,  rad/sec 

Components  of  S^ss  about  X-,  Y-,  and  Z-body  axes,  respectively,  rad/sec 
Dynamic  pressure,  pV2/2,  lb/ft2 


59 


AEDC-TR-77-1  26 


Reg  Reynolds  number  based  on  model  length 

RDELCG  Moment  arm  for  X axis  center-of-gravity  transfer,  ft 

rpm  Revolutions  per  minute 

S Wing  reference  area,  ft2 

T Engine  thrust,  lb 

Tx  Component  of  thrust  along  X-body  axis,  lb 

Ty  Component  of  thrust  along  Y-body  axis,  lb 

Tz  Component  of  thrust  along  Z-body  axis,  lb 

t Time,  sec 

u,v,w  Components  of  total  velocity  along  X-,  Y-,  and  Z-body  axes,  respectively, 

ft/sec 

V Total  velocity,  ft/sec 

X,Y,Z  Body  axes 

x,y,z  Linear  distance  along  X-,  Y-,  and  Z-body  axes,  respectively,  ft 

a Angle  of  attack,  deg 

(3  Angle  of  sideslip,  deg 

7 Flight  path  angle,  deg 

5 Angle  between  rotation  vector  ^ and  the  velocity  vector,  deg 

8a  Aileron  deflection,  positive  when  trailing  edge  of  right  aileron  is  down,  deg 

8h  Elevator  deflection,  positive  when  trailing  edge  is  down,  deg 

6r  Rudder  deflection,  positive  when  trailing  edge  is  left,  deg 

e Angle  between  thrust  line  and  X-body  axis,  deg 

6 Angle  between  X-body  axis  and  horizontal  measured  in  vertical  plane,  deg 
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p Air  density,  slug/ft3 

r Time  constant 

At  Incremental  time 

0 Angle  between  Y-body  axis  and  horizontal  measured  in  vertical  plane,  deg 

0 e Engine  rotor  angular  velocity,  rad/sec 

0 Angle  between  Y-body  axis  and  vertical  measured  in  horizontal  plane,  deg 

^ Resultant  angular  vector,  rad/sec 

ftb/2V  Nondimensional  rotation  rate 

tls&  Steady-state  angular  velocity  (component  of  Tl  along  velocity  vector),  rad/sec 

SUBSCRIPTS 

rb  Indicates  rotation-balance  data 

sp  Short  period 

SUPERSCRIPT 

Derivative  with  respect  to  time 
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